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KR 553431 Bl KO Previously issued patent KR 2004091271 

Image segmentation optimum threshold value finding method for use in 
image recognition, involves gaining minimum entropy value related to 
gray level as threshold value by using fixed point iteration 
based on entropy values 

Original Titles: 

. . .Method and apparatus for finding optimal threshold for image 
segmentation 

. ..Procede et dispositif d'etablissement d'un seuil optimal pour la 
segmentation d ' image 

. . .Method for finding optimal threshold for image segmentation 
Inventor : SHIN Y 

. . . SHIN Y s 

Alerting Abstract . . .NOVELTY - The method involves computing entropy 
values corresponding to gray values in a histogram . Entropy values of 
gray levels related to possible optimal thresholds are obtained. A 
threshold value is obtained by comparing entropy values and selecting a 
minimum entropy value. The minimum entropy value related to the 
gray level is gained as the threshold value by using a fixed point 
iteration based on the computed entropy values. USE - Used for finding an 
optimum threshold value for an image segmentation in an image 
recognition. . . 

. . .ADVANTAGE - The method analyzes an entropy characteristic of image based 
on fixed point iteration , thereby effectively and quickly finding an 
optimal threshold for image segmentation of an image . 



. . .DESCRIPTION OF DRAWINGS - The drawing shows a flowchart for explaining 
the step of gaining a gray level corresponding to the minimum 
entropy by using fixed point iteration . 

Title Terms ... /Index Terms /Additional Words: THRESHOLD / 
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Shin, Yong-Shik . . . 

. . . SHIN Y S . . . 



2 



10817551 Method for Finding Optimal Threshold for Image Segmentation 



. . . shin y s . . . 

. . . Shin, Yong-Shik 
Original Abstracts: 

A method and apparatus for finding the optimal threshold for image 
segmentation in image recognition is disclosed. The method includes the 
steps of: a) gaining histogram distribution of an image; b) computing 
entropy values corresponding to gray levels in the histogram ; and c) 
gaining a minimum entropy value corresponding to the gray level as 
the threshold value by using a fixed point iteration FPI based on 
the computed entropy values... 

. . .A method for finding the optimal threshold for image segmentation 
in image recognition is disclosed. The method includes the steps of: a) 
gaining histogram distribution of an image; b) computing entropy values 
corresponding to gray levels in the histogram ; and c) gaining a 
minimum entropy value corresponding to the gray level as the 
threshold value by using a fixed point iteration FPI based on the 
computed entropy values. 
Claims : 

A method for finding a threshold value in image segmentation , the 
method comprising the steps of: a) gaining histogram distribution of an 
image; b) computing entropy values corresponding to gray levels in the 
histogram ; and c) gaining a minimum entropy value corresponding to the 

gray level, as the threshold value by using a fixed point 
iteration FPI based on the computed entropy values... 

...What is claimed is: <b>l</b>. A method for finding a threshold value 
in image segmentation , the method comprising the steps of: a) gaining 
histogram distribution of an image; b) computing entropy values 
corresponding to gray levels in the histogram ; and c) gaining a 
minimum entropy value corresponding to the . gray level as the 
threshold value by using a fixed point iteration FPI based on the 
computed entropy values. 
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Original Publication Data by Authority 



Original Abstracts : 

. . .of the lung fields and costophrenic angle (CP) regions in 
posteroanterior (PA) chest radiographs, wherein image segmentation 
based on gray - level threshold analysis ( S3 , 1003 ) is performed by 
applying an iterative global gray - level thresholding method ( S5 , 
1005 ) to a chest image based on the features of a global gray - level 
histogram (S3, 1003 ) . Features of the regions in a binary image 
constructed at each iteration are identified and analyzed to exclude 
regions external to the lung fields. The initial lung contours that result 
from this global process are used to facilitate a local gray - level 
thresholding method ( S6 , 1006 ) . Individual regions-of -interest (ROIs) 
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are placed along the initial contour. A procedure is implemented to 
determine the gray - level thresholds to be applied to the pixels 
within the individual ROIs. The result is a binary image, from which... 

. . .of the lung fields and costophrenic angle (CP) regions in 
posteroanterior (PA) chest radiographs, wherein image segmentation 
based on gray - level threshold analysis is performed by applying an 
iterative global gray - level thresholding method to a chest image 
based on the features of a global gray - level histogram . Features 
of the regions in a binary image constructed at each iteration are 
identified and analyzed to exclude regions external to the lung fields... 
...initial lung contours that result from this global process are used to 
facilitate a local gray - level thresholding method. Individual 
regions-of -interest (ROIs) are placed along the initial contour. A 
procedure is implemented to determine the gray - level thresholds to be 
applied to the pixels within the individual ROIs. The result is a 
binary image, from which final contours are constructed. Smoothing... 

...is employed on a column -by- column basis to identify initial diaphragm 
points, and maximum gray - level information is used on a row-by- row 
basis to identify initial costal points. Analysis of initial diaphragm and 
costal points allows for appropriate... 

. . .automated segmentation of the lung fields and costophrenic angle (CP) 
regions in posteroanterior (PA) chest radiographs , wherein image 
segmentation based on gray - level threshold analysis is performed by 
applying an iterative global gray - level thresholding method to a 
chest image based on the features of a global gray - level histogram . 
Features of the regions in a binary image constructed at each iteration 
are identified and analyzed to exclude regions external to the lung 
fields. The initial lung contours that result from this global process are 
used to facilitate a local gray - level thresholding method. Individual 
regions-of -interest (ROIs) are placed along the initial contour . A 
procedure is implemented to determine the gray - level thresholds to be 
applied to the pixels within the individual ROIs. The result is a 
binary image, from which final contours are constructed. Smoothing 
processes are applied, including a... 

...is employed on a column-by-column basis to identify initial diaphragm 
points, and maximum gray - level information is used on a row-by- row 
basis to identify initial costal points . Analysis of initial diaphragm 
and costal points allows for appropriate adjustment of CP angle ROI . . . 

...segmentation of the lung fields and costophrenic angle (CP) regions in 
posteroanterior (PA) chest radiographs wherein image segmentation 
based on gray - level threshold analysis is performed by applying an 
iterative global gray - level thresholding method to a chest image 
based on the features of a global gray - level histogram . Features of 
the regions in a binary image constructed at each iteration are 
identified and analyzed to exclude regions external to the lung fields. 
The initial lung contours that result from this global process are used to 
facilitate a local gray - level thresholding method. Individual 
regions-of -interest (ROIs) are placed along the initial contour. A 
procedure is implemented to determine the gray - level thresholds to be 
applied to the pixels within the individual ROIs. The result is a 
binary image, from which final contours are constructed. Smoothing 
processes are applied, including a unique... 
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...is employed on a column-by-column basis to identify initial diaphragm 
points, and maximum gray - level information is used on a row-by-row 
basis to identify initial costal points . Analysis of initial diaphragm 
and costal points allows for appropriate adjustment of CP angle ROI 
positioning. . . 

...of the lung fields and costophrenic angle (CP) regions in 
posteroanterior (PA) chest radiographs, wherein image segmentation 
based on gray - level threshold analysis is performed by applying an 
iterative global gray - level thresholding method to a chest image 
based on the features of a global gray - level histogram . Features 
of the regions in a binary image constructed at each iteration are 
identified and analyzed to exclude regions external to the lung 
fields. The initial lung contours that result from this global process are 
used to facilitate a local gray - level thresholding method. Individual 
regions-of -interest (ROIs) are placed along the initial contour. A 
procedure is implemented to determine the gray - level thresholds to 
be applied to the pixels within the individual ROIs. The result is a binary 
image, from which final contours are constructed. Smoothing processes 
are applied, including a unique adaptation of a rolling. . . 

...is employed on a column-by-column basis to identify initial diaphragm 
points, and maximum gray - level information is used on a row-by- row 
basis to identify initial costal points. Analysis of initial diaphragm 
and costal points allows for appropriate adjustment of CP angle ROI 
positioning. Polynomial curve- fitting is... 

. . .of the lung fields and costophrenic angle (CP) regions in 
posteroanterior (PA) chest radiographs, wherein image segmentation 
based on gray - level threshold analysis is performed by applying an 
iterative global gray - level thresholding method to a chest image 
based on the features of a global gray - level histogram . 
Features of the regions in a binary image constructed at each 
iteration are identified and analyzed to exclude regions external to the 

lung fields . The initial lung contours that result from this global 
process are used to facilitate a local gray - level thresholding method. 
Individual regions-of -interest (ROIs) are placed along the initial contour. 
A procedure is implemented to determine the gray - level thresholds to 

be applied to the pixels within the individual ROIs. The result is a 
binary image, from which final contours are constructed . Smoothing 
processes are applied, including a unique adaptation of a rolling ball 
method. CP angles are... 

...is employed on a column-by-column basis to identify initial diaphragm 
points, and maximum gray - level information is used on a row-by-row 
basis to identify initial costal points. Analysis of initial diaphragm and 
costal points allows for appropriate adjustment of CP angle ROI 
positioning. Polynomial curve-fitting is used to combine the diaphragm... 

. . .of the lung fields and costophrenic angle (CP) regions in 
posteroanterior (PA) chest radiographs, wherein image segmentation 
based on gray - level threshold analysis (S3, 1003) is. performed by 
applying an iterative global gray - level thresholding method (S5, 1005) 
to a chest image based on the features of a global gray - level 
histogram (S3, 1003). Features of the regions in a binary image 
constructed at each iteration are identified and analyzed to exclude 
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regions external to the lung fields. The initial lung contours that 
result from this global process are used to facilitate a local gray - 
level thresholding method ( S6, 1006) . Individual regions -of -interest 
(ROIs) are placed along the initial contour .A procedure is implemented to 
determine the gray - level thresholds to be applied to the pixels 
within the individual ROIs. The result is a binary image, from which final 
contours are constructed. 



...champs pulmonaires et des angles costo-phreniques, dans des 
radiographics postero-anterieures des poumons . Cette segmentation d' 
image s'effectue d'apres une analyse de seuil des niveaux de gris (S3, 
1003) , par. . . 

...d'une technique de seuillage global et iteratif (S5, 1005) des niveaux 
de gris (5, 1005 ), sur une image des poumons, d'apres les 
caracteristiques d'un histogramme global (S3, 1003) des niveaux de 
Claims : 

...of the chest region; andconstructing, based on said lung fields 
identified in said processed image , first initial lung segmentation 
contours for said posteroanterior chest image. 



. . .of the chest region; andconstructing, based on said lung : fields 
identified in said processed image , first initial lung segmentation 
contours for said posteroanterior chest image... 

lung segmentation contours for said posteroanterior chest 

image; per forming, based on said first initial lung segmentation contours 
, local threshold analysis to construct second initial lung segmentation 
contours for said posteroanterior chest image; andapplying a rolling ball 
filter to said second initial lung segmentation contours to smooth the 
shape of said second initial lung segmentation contours. 
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Method of automatically segmenting lung radiographs 
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Alerting Abstract ...and posterior margins are determined (1003,1004) . A 
Sobel filter (1005) is applied and global gray - level thresholding is 
performed (1006). The resulting contour is smoothed (1007), e.g. using 
running mean and rolling ball techniques. Adaptive local gray - level 
thresholding (1009) is followed by further smoothing (1010) and third 
order, least-squares polynomials are... 
...1006 Global gray - level thresholding... 

...1009 Adaptive local gray - level thresholding... 

Original Publication Data by Authority 



Original Abstracts : 

...for the automated segmentation of the lung regions in lateral chest 
radiographs (10)- based on gray - level threshold analysis. Approximate 
outer bounds on the extent of the lung fields in the image are identified 
to restrict the region further analyzed (16) . An iterative global gray - 
level thresholding method ( 20 ) is applied based on the features of a 
global gray - level histogram . Features of the regions in a 
binary image constructed at each iteration are identified and subjected to 
a modified. . . 

...lung field. Individual regions-of -interest (ROIs) are placed along the 
initial contour. The single gray - level threshold to be applied to 
the pixels within the individual ROIs is determined (1009) . A final contour 
is constructed. . . 

. . .This is achieved according to the invention by providing an improved 
computerized, automated method for image segmentation based on gray - 
level threshold analysis. A unique method for identifying an 
approximate outer bounds on the extent of the lung fields in the image is 
performed to restrict the region further analyzed. An iterative global 
gray - level thresholding method is applied based on the features of a 
global gray - level histogram . Features of the regions in a 
binary image constructed at each iteration are identified and subjected ' to 

a modified analysis to exclude regions external to the lung field. The 
initial lung region contour that results from this global process is used 
to facilitate a novel adaptive local gray level thresholding method. 
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Individual regions - of -interest (ROIs) are placed along the initial 
contour. The dimensions of the several ROIs are... 

. . .upon the patient anatomy enclosed therein. A unique procedure is 
implemented to determine the single gray - level threshold to be applied 

to the pixels within the individual ROIs. A composite binary image 
results, and a final contour... 

...for the automated segmentation of the lung regions in lateral chest 
radiographs (10) based on gray - level threshold analysis. Approximate 
outer bounds on the extent . of the lung fields in the image are 
identified to restrict the region further analyzed (16) . An iterative 
global gray - level thresholding method (20) is applied based on the 
features of a global gray - level histogram . Features of the regions 
in a binary image constructed at each iteration are identified and 
subjected to a modified analysis to exclude regions external to the 
lung field. Individual regions-of -interest (ROIs) are placed along the 
initial contour. The single gray - level threshold to be applied to the 
pixels within the individual ROIs is determined (1009) . A final contour 
is constructed to enclose "on" regions (26... 

...chaque iteration les caracteristiques des regions dans une image binaire 
construite, puis on les soumet a une analyse modifiee de facon a 
exclure les regions externes au champ du poumon. On place le . . . 
Claims : 

. . .posterior margins in said second image data to produce third image 
data performing iterative global gray - level thresholding on said third 
image data to identify a first initial lung segmentation contour; 
and smoothing said first initial lung segmentation contour to produce a 
second initial lung segmentation contour. 
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Method and apparatus for finding optimal threshold for image 
segmentation 

Procede et dispositif d 1 etablissement d'un seuil optimal pour la 

segmentation d 1 image 
INVENTOR : 

Shin, Yong-Shik . . . 

. . .ABSTRACT A3 

A method and apparatus for finding the optimal threshold for image 
segmentation in image .recognition is disclosed. The method includes 
the steps of: a) gaining histogram distribution of an image; b) 
computing entropy values corresponding to gray levels in the 
histogram ; and c) gaining a minimum entropy value corresponding to 
the gray level as the threshold value by using a fixed point 
iteration FPI based on the computed entropy values. 

...SPECIFICATION of the Invention 

The present invention relates to a method and apparatus for finding a 
threshold for image segmentation ; and, more particularly, to a 
method and apparatus for finding the optimal threshold for image 
segmentation in image recognition. 

Description of Related Arts 

Generally, a process of finding the optimal threshold for the image 
segmentation is fundamental and important process in the image 
recognition. The process of recognition is necessary to distinguish an 
object from a background of an image. 

The optimal threshold can be found easily based on a bimodal type 
histogram distribution graph and in above case, it is located at a 
lowest point of histogram distribution curve. There are many methods 
introduced for finding the optimal threshold . 

A first method is stochastic method to find the optimal threshold . 
That is, a histogram distribution of an image is assumed as the bimodal 
type and a gray level having the minimum sum of variance is select as 
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the optimal threshold . A second method finds the optimal threshold 
based on the Shannon entropy. A gray level having the minimum 
entropy is selected as the optimal threshold . There is also another 
method using fuzziness during computing entropy of image. This method 
selects a gray level having minimum fuzziness as the optimal 
threshold . 

Fig. 1 is a flowchart explaining a conventional method for finding the 
optimal threshold . Referring to Fig. 1, at step S101, a histogram 
distribution of an image is computed. Entropies of all gray levels 
distributed in the histogram are computed at step S102. All computed 
entropies are compared one another and a gray level having the 
lowest entropy is selected at step 103. However, a process time of the 
conventional method is increased. . . 

...an object of the present invention to provide a method for effectively 
finding an optimal threshold for image segmentation of an image 
having multi thresholds by analyzing entropy characteristic of image 
based on a fixed point iteration and fuzzy entropy . 

In accordance with an aspect of the present invention, there is 
provided a method for finding a threshold value in image 
segmentation , the method including the steps of: a) gaining histogram 
distribution of an image; b) computing entropy values corresponding to 
gray levels in the histogram / and c) gaining a minimum entropy 
value corresponding to the gray level as the threshold value by 
using a fixed point iteration FPI based on the computed entropy 
values . 

Brief Description of the Drawing (s) 
The above and. . . 

. . . in which: 

Fig. 1 is a flowchart explaining a conventional method for finding the 
optimal threshold ; 

Fig. 2 is a flowchart for explaining a method for finding the optimal 
threshold for image segmentation in accordance with a preferred 
embodiment of the present invention; 

Fig. 3 is a graph. . . 
. . .accordance with a preferred embodiment of the present invention; 

Fig. 4 is a graph showing gray level distribution curve for 
obtaining a gray level of minimum entropy in accordance with a 
preferred embodiment of the present invention; 

Fig. 5 is a flowchart for explaining the step S203 in Fig. 2 for ■ 
gaining a gray level corresponding to the minimum entropy by 
using FPI in accordance with a- preferred embodiment of the present 
invention; and 

Fig. 6 is a. . . 

. . .forth hereinafter. 

Fig. 2 is a flowchart for explaining a method for finding the optimal 
threshold for image segmentation in accordance with a preferred 
embodiment of the present invention. 

Referring to Fig. 2, a histogram distribution of an image is obtained 
at step 201. Entropies of gray levels are computed at step S202. 
After computing entropy value at step S202, a gray level of minimum 

entropy is gained by using the fixed point iteration ( FPI ) 
based on the computed entropy value at step S203 . 

In the step S202, the entropy of the gray level is computed by 
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measuring fuzzy entropy of corresponding gray level . Hereinafter, 
the computation of fuzzy entropy is explained in detail. 

If there is an M x N size of image I having L gray levels , a gray 
level of pixel (x,y) is defined as I(x,y) and (mu) I) ) ( Ii , j ) ) ) 
represents fuzziness of gray scale of pixel (x,y) . Therefore, the 
image I can be expressed as following equation. , wherein. . . 

...j))) ^ l;i = 0,1,..., M-l; j = 0,1,..., N-l 

If a gray level g has a frequency of generation h(g) in entire 
image I then an average gray level (mu)0)) of a background can be 
expressed as following equation 2 and an average gray level (mu)D) 
of an object can be expressed as following equation 3. 

The averaige gray levels (mu)D) and (mu)o)) can be considered as 
target values of the threshold value T. That is, the fuzziness can be 
expressed as a difference between gray level of a pixel (x,y) and a 
gray level of a region including the corresponding gray level . 
Therefore, the difference of gray levels is smaller, as larger as the 
fuzziness is. Gray levels of all pixels in an image for a threshold 
T must have certain fuzziness either to an object or background. The 
fuzziness of a pixel can be expressed as following equation. 

In a meantime, when a gray level of a certain pixel is included in 
a specific region, the fuzziness must to be... 

. . .as a fuzziness of fuzzy set and there are various entropy-equations 

disclosed for computing, fuzzy entropy . If the entropy equation of one 
independent variable is expanded to 2 dimensional image region. . . 

...decreased in a region of (0.5, 1) . In a meantime, if fuzziness of all 
gray level included in the image are about 0.5 then entropy E(I) has 1 
as the maximum value. 

After obtaining entropy values at step S202, a gray level of 
minimum entropy is obtained by using a fixed point iteration ( FPI 
) method at step S203. 

Fig. 4 is a graph showing gray level distribution curve for 
obtaining a gray level of minimum entropy in accordance with a 
preferred embodiment of the present invention. 

Referring to Fig. 4, obtaining a gray level of minimum entropy 
is explained hereinafter. 

Fig. 5 is a flowchart for explaining the step S203 in Fig. 2 for 
gaining a gray level corresponding to the minimum entropy by 
using FPI in accordance with a preferred embodiment of the present 
invention. 

At step 501, possible optimal threshold values Pi)) for obtaining a 
gray level of minimum entropy are obtained based on the graph of 
Fig. 4. 

After obtaining the Pi) ) , gray levels are sequentially obtained 
from left to right Pi) ) . 

At step 503, an optimal threshold of gray level having minimum 
entropy is obtained by comparing entropy values of gray levels of 

Pi) ) . 

Fig. 6 is a flowchart for explaining step 501 of Fig. 5 in... 

. . .Gmin) ) , Gmax) ) , Pi and gcal) ) are set as follows, gmin) ) is set as 
possible minimum gray level by selecting a lowest value of a gray 
level distribution curve on Fig. 4, and gmax)) is set as possible 
maximum gray level by selecting a highest value of a gray level 
distribution curve on Fig. 4. And Gmax)) is set as equal to gmax)) and 
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Gmin. . . 

...steps 602 and 608 are reputedly performed. For helping to understand 
steps for obtaining optimal threshold of Fig. 6, pseudo code is shown 
in below table . 

As mentioned above, the present invention can quickly find the optimal 
threshold value by analyzing entropy characteristic of image based on 
a segmentation completion condition and a fixed point iteration . 

While the present invention has been described with respect to certain 
preferred embodiments, it will... 

. . . CLAIMS A2 

1. A method for finding a threshold value in image segmentation , 

the method comprising the steps of: 

a) gaining histogram distribution of an image; 

b) computing entropy values corresponding to gray levels in the 
histogram / and 

c) gaining a minimum entropy value corresponding to the gray 
level as the threshold value by using a fixed point iteration 

FPI based on the computed entropy values. 

2. A method as recited in claim 1, wherein... 

. . .c-1) obtaining a plurality of possible optimal thresholds; 

c-2) obtaining entropy values of gray levels corresponding to the 

obtained possible optimal thresholds; and 
c-3) obtaining the threshold value by comparing entropy values and 

selecting minimum entropy value. 

3. A method as recited in claim 2, wherein each of the possible optimal 

thresholds is obtained by obtaining a value of possible maximum gray 

level having maximum entropy value, a value of possible minimum 
gray level having minimum entropy value and obtaining possible 
optimal threshold by adding two values of the possible maximum 
gray level and the possible minimum gray level and dividing 
the sum of addition by half. 

4. A method as recited in claim... 

. . .possible optimal thresholds are obtained by changing one of the value of 
the possible maximum gray level and the value of the possible 
minimum gray level according to comparison of entropy values of 
the possible maximum gray level , the possible minimum gray 
level and obtained optimal threshold and by newly obtaining a 
possible optimal threshold based on the changed values of the 
possible maximum gray level and the value of the possible minimum 
gray level . 

5. A method as recited in claim 2. wherein the step c-1) includes the 

steps of: 

c-i) obtaining an initial possible optimal threshold , an initial 

possible maximum gray level having maximum entropy value and an 
initial possible minimum gray level having minimum entropy 
value by setting Gmin) ) to have the initial possible minimum gray 
level , setting Gmax) ) to have the initial possible maximum gray 
level , setting gmin) ) and gmax) ) to have identical values Gmin) ) and 
Gmax) ) , respectively for not influencing change of value of Gmin) ) 
and Gmax) ) , setting Pi) ) to have the initial possible optima 
threshold by computing equation Pi) ) =( (gmin) ) + gmax) ) ) /2) and 
setting gcal) ) to have the identical value... 
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.. .value of gf ix) ) to have the value of Gmin) ) / 

c-vi) obtaining new possible optimal threshold pi)) based on changed 
value of gmin)) and gmax) ) by an equation as: Pi)) = (gfix. . . 

. . . s ; 

c-ix) if there are not identical two Pi))s, determining next possible 
optimal threshold by setting gtemp) ) to have the value of pi+1) ) 
and setting gcal) ) to have. . . 

. . .c-viii) ; and 

c-x) if there are identical any two Pi))s, selects the threshold value 
by comparing entropy values of corresponding Pi))s and selecting Pi 
having minimum entropy value as the threshold value . 
6. Apparatus for finding a threshold value in image segmentation , 
the said apparatus comprising: 

(a) means for gaining a histogram distribution of an image; 

(b) means for computing entropy values corresponding to grey levels 
in a histogram / and 

(c) means for gaining a minimum entropy value corresponding to a 
grey level as a threshold value; whereby the said minimum 
entropy value is gained by using a fixed point interaction based on 
computed entropy values . 
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signals at the sub-subsets) obtained for any split is less than a 
selected threshold , the subset is not split. Moreover, if the number of 
subsets reaches a selected limit... 



..maximum benefit is obtained with a fixed number of subsets, the subset 
with the highest entropy is split in each iteration . 

In the method described thus far, the candidate questions were 
limited to those of the form "Is the... 
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